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Abstract 

The solution of the Muskhelishvili-Omnes Integral Equation is ambiguous by 
a real polynomial. The coefficients of this polynomial can be fixed either by the 
knowledge of the low energy parameters or by the asymptotic behavior of the form 
factor. The role of the contact terms of the low energy effective Lagrangian is 
explicitly analysed. 
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There have been discussions of contact terms in the phenomenological models de- 
scribing low energy pion phenomena 0, |3[]. They seem to play an important role in 
analysing the pion form factor and the process 771"° — > 7r + 7r~. The main difference between 
this model and the usual vector meson model is the existence of a contact term where 
the photon interacts directly with the charged pions and not by the intermediate of the 
vector p meson. Phenomenological analysis of this contact term in the form factor has 
recently been given 0. 

The purpose of this note is to analyze the existence of the contact term within the 
spirit of the more general approach to the form factor problem using the Muskhelishvili- 
Omnes (MO) integral equation ||. It is based on the assumption of the elastic unitarity 
and dispersion relation for the pion form factor. The main conclusion of our analysis is 
that the existence of the contact term is in conflict with unitarity in the form of the phase 
theorem. Its existence can, however, be simulated by a polynomial ambiguity contribution 
to the solution of the MO equation. This polynomial is interpreted physically as either the 
correction due to the inelastic contribution to the unitarity relation or the modification 
of the asymptotic behavior of the form factor. 

Let us first discuss the solution of this problem from a more general viewpoint. The 
form factor V(s) is an analytic function with a cut from 4m^ to oo. 

Assuming that V(s) is polynomially bounded and that V(s)s~( n+1 ' — > as s — > oo, 
n > 0, we can write the following subtracted dispersion relation: 

, _ s n+1 r°° ImV(z)dz 
V(s) = l + a lS + ...a n s n + / 1 

71 J Ami (Z) n+i {Z — S — 16) 

where we have used the Ward identity V(0) = 1. Around s — 0, the dispersion integral 
is of the order (s) n+1 and can be neglected, hence the low energy theorem is recovered. 
Needless to say, a n are, apart from a factorial factor n\, the derivatives of V(s) evaluated 
at s = 0. These could be taken from experiments or could be determined by sum rules 
using some assumptions on the asymptotic behavior of the form factor. 

For example, if we assume that as s — > oo, V(s)/s — > 0, a once subtraction dispersion 
relation for V(s) can be written; the first and higher derivatives of V(s) can be written 
as: 

n! f°° ImV(z)dz 

7T J Ami Z*' 

The mathematical problem is now clear: Find the solution of the integral equation of 
the MO type for V(s) with its imaginary part given by the elastic unitarity: 

ImV(s) = V{s)e- i5{s) sin 5{s) (3) 

with the boundary conditions around s = given by Eq. (|1|) and where 5 is the the strong 
P-wave I = 1 7T7T phase shifts. 

To solve this integral equation Eq. ([I]) which is of MO type M, let us define the 
function O(s) normalized to unity at s = 0: 

r^/ \ / S f°° S(z)dz . , A . 

n (s) = exp - / — 4 

7T J Ami Z(Z — S — ie) 



v""»(o) = n 4 r. -^r <2) 
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The solution for our integral equation is: 

V(s) = P n (s)Q(s) (5) 

where P n (s) is a polynomial in s of order n with real coefficients, P„(0) = 1. From eq. 
(ID, the phase of the form factor V(s) is 5(s) which is also a direct consequence of the 
unitarity relation ||. It will be called in the following as the phase theorem. 

The ambiguity of the solution of the MO equation due to the polyomial is not surprising 
because its solution depends on the knowledge of the asymptotic behavior of the form 
factor. The polynomial could be interpreted as extra parameters which are introduced in 
the problem without violating the phase theorem; in many cases they can be determined 
using experimental data at low energy. 

This is indeed true by expanding the function Q(s) in a power series in s, and compare 
it with the derivatives given in Eq. ([5]) and Eq. ([J). The expansion of Q(s) in the Taylor's 
series around s = is possible because fl(s) is an analytic function with a cut from 4m^ 
to oo. 

In the special case where only two terms in the series are known, i.e. the pion charge 
by Ward's indentity and its first derivative a% which is proportional to the r.m.s. pion 
radius by experiment, the solution of our integral equation is given by: 

V(s) = {1 + s(a 1 -fi'(0))}fi(s) (6) 
where fl (0) denotes the first derivative of f2(s) evaluated at s = and is given by: 

= i J- ^ 

7T J Ami z 

The presence of the term fi'(0) is to ensure the boundary condition for V(s) is satisfied. If 
n (0) = there is no need to introduce a momnomial term in the form factor and hence 
the asymptotic behavior of the form factor is solely determined by Q(s). (Using partial 
fraction for the dispersion integral together with the sum rules Eq. (y) the asymptotic 
behavior of V(s) can be determined). 

It is straightforward to generalise the solution of Eq. ([!]) for other values of n. For 
example when n = 2, the solution for the integral equation Eq.( [1]) is obtained by adding 
to the curly bracket on the righthand side of Eq. (§) a term: 

s 2 (a 2 - ai n'(0) - a ^^ + fi' 2 (0)) (8) 

Let us now look at the experimental data. For a given set of P-wave hit phase shifts 
5, one can, in principle, construct the function Q(s) as given by Eq. (||) by the numerical 
method. We prefer instead to construct this function using the inverse amplitude method 
or the [0,1] Pade approximant method for the one loop Chiral Perturbation Theory for the 
pion form factor and verify that the resulting phases are in agreement with experiments 
H]. The function Q(s) can be parametrised as follows: 

fi(s) = 1 - s/s R - 7^{(s - K)^(s) + 2s/3} (9) 



2 



where /„- = 0.093GeV, and sr is related to the p mass squared m 2 by requiring that the 



real part of the denominator of Eq. @ vanishes at the p mass; H„ 
integral over the phase space factor: 



:s) is a well-known 



(2 - 2^1 - Aml/s) lot 



Anil 



2m, 



+ in J 1 — 4m 2 /s 



(10) 



for s > 4m 2 ; for other values of s, H nn (s) can be obtained by analytic continuation. Q(s) 
as given by Eq. @ has a ghost at a very large negative value of s with a small residue 
and is irrelevant for our low energy calculation. 

Although we do not try the best fit to the experimental data on the P-wave ixix phase 
shifts, an adequate agreement is obtained between theory and experiment as shown in 
Fig. (1) where m p = 0.775GeV. 

Eq. @ can also be obtained by the vector meson dominace model with the strong 



pixix coupling constant given by the KSRF relation [10 



If it was assumed that there is no polynomial ambiguity P n (s) = 1, we would have: 



V(s) = n(s) 



11) 



Fitting this expression with the physical value of the p mass, the pion rms radius 
calculated by this expression is too low by more than 10% and the maximum value of 
the form factor squared at the p peak (the p leptonic width) is too low by 30%, Fig. (2). 
This implies that the assumed elastic unitarity relation is not good enough. One can 
either rewrite the integral equation taking into account of the inelastic effect which was 



previously studied [DJ or to make more subtractions in the integral equation and keeping 
the elastic unitarity relation, in order to minimize the high energy contribution. We choose 
first the latter method. (The two-loop inverse amplitude or the [0,2] Pade approximant 
method are given recently and are found to give excellent agreement with the pion form 
factor data [§]). 

Let us write a twice subtracted dispersion relation, including the inputs as the pion 
charge and the experimental value of the pion radius. The integral equation now reads: 



V(s) = 1 + ais H 

7T J 4ml 



V(z)e~ i5 ^ sin 5(z)dz 



z z (z — s — le) 



(12) 



where a\ = l/6(< r 2 >) with < r 2 > is the rms squared radius of the pion. The factor z 2 
in the denominator of the dispersion integral suppresses the high energy contribution. 

. the solution of the 



Using the experimental value < r >= 0.439 ± .008/m z 
integral solution is given by Eq. (Bf) and is numerically equal to |L2L |22 



V(s) 



0.15s /m 2 p )Q(s 



(13) 



The p leptonic width and the square of the modulus of the pion form factor are now in 
agreement with the experimental data, Fig. (2). 
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The other possibility to get agreement with experimental data is to introduce the 
inelastic effect which could occur at higher energy. The integral equation now reads: 



v(s) = 1 + ir dz rmm^> (14) 

7T J Ami Z \ Z — S — It) 

where f(s) = (2i)~ 1 (r]e 2i<5 ^- ) — 1), rj is the inelastic factor (rj = 1 for elastic tctc scattering) 
and a is the inelastic spectral function due to the contribution of higher intermediate 
states in the unitarity relation. The solution of this integral equation is known: 

■t t t \ , , M s [°° , 2Re(a(z)e i5 &) 1 

v{.) = n(.){i + - 1 dz (l + vi ' z) \J_ tm n( z) z{z - s - u) ] (15) 

Si is the inelastic threshold and for the pion form factor problem one can take it to be 
around 1 — 1.3GeV 2 . In reality, it could be as high as 2 GeV 2 in the p' resonance where 
the inelastic effect would become important. 

The dispersion integral in Eq. (15) is seen to be an analytic function with a cut 



starting at s, to oo and hence it can be expressed as a Taylor's series around s = 0. The 
radius of convergence of this series is Sj. Hence the RHS of Eq. (|i~5|) is simply the solution 
of the MO equation as given by Eq. (Bl) and the fit to the experimental data is given by 



Eq. (|T1). 

For a feel for the inelastic effect, let us see how does the p'(1500) , which dominantly 
decays into 4tt, influence the pion form factor. The inelastic spectral function can then be 
approximated by a delta function in the dispersion integral which yields a term propor- 
tional to Q(s) multiplying with the p' propagator. This assures the validity of the phase 
theorem below the inelastic threshold and that neither the p' propagator alone nor a pure 
contact term can be written. For s in the p region, one can expand the p propagator in 
a series in s and fit to the experimental data to get the desired result Eq. (|13|). 

Let us now compare our result with models with a contact term 0] and, in particular, 
the hidden local symmetry model ||. At the tree level, the hidden local symmetry model 
gives: 



= "5 + 1 + I^S < 16 > 

which can be rewritten as: 

utq, , m \ + s(a/2 - 1) , s 

V HLS {s) = -£ ^ (17) 



In order to make the p unstable we could try to replace the factor m^/lm 2 , — s) in Eq. 
(|16|) by fl(s) leaving the contact term unchanged. This would violate the phase theorem 
unless the constant term cancels each other to give a net result which is equal to fi(s); 
this solution does not fit with the experimental data as explained above. 

We should use instead Eq. fllTD and interprete it as the limit large s, V HLS — > (1 — a/2) 
(in the limit of zero width). Then the hidden symmetry solution for the vector form factor 
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must contain a polynomial Pi(s) or: 



yHLS = {l + ZUL_ S ^ S 



a/2 - 1 



m 2 



Compare this equation with Eq. (|13|) we have a = 2.30. Within the validity of the above 
interpretation of the contact term, this result is in agreement with that obtained by the 
authors of the reference ||] who begin with Eq. (|16|), then make the p unstable. The 
unitarisation scheme is done with the N/D method for the P-wave irn scattering. 

The limit of large value of s of V HLS should be taken with some reservation, because 
the inelastic effect interpretation of the polynomial ambiguity is only valid for s < Sj, Eq. 



(15|) and Sj is experimentally not very large. 

Another example of the contact term [p], [2], |J is discussed in the jtt — > 7r + 7r~ process, 
in particular, the hidden symmetry model ||. In the VMD model, the 77?° — > 7r + 7r~ is 
written as: 

A*"{s, t, u) = + + (19) 

3 m l p — s m A p — t m l p — u 

where s, t, u are the standard invariant kinematics. With a contact term it can be written 
as: 

, A rn 2 n m\ m 2 n 

A 3 *(s, t, u) = ^{-^ + - V4 + ~ c > (2°) 

3 — c m p — s m p — t m p — u 

where c is proportional to the strength of the contact term and A is given by the anomaly: 

A = -4tt = 9.7GW- 3 (21) 



47T2/3 

Eq. (]20D can be rearranged to give: 



\,r ml c s 



A**(a, t, u) = £{[-5-Ml + ^—^)] + [s-t} + [s^ u]} (22) 



3 LL m p — s 3 — cm p ' 



Eq. ( p!9|) yields a decay width T(p — > 7r7) = 36KeV which is too small compared with 
the experimental data. Eq. ( ^0|) with c = 1 yield, on the other hand, a width which is 



9/4 times larger and appears to be in better agreement with the data [19 



To make the p unstable, we could make a substitution m 2 p /(m 2 p — s) with Q(s) in 
Eq. ([B]) and Eq. fl2"0p, then both resulting equations would violate the phase theorem. 
Following the spirit of the form factor calculation, one could try to make this replacement 
in Eq. (|2~2D but the resulting equation still violates the phase theorem when the P-wave 
amplitude is projected out from the full amplitude. How to get a correct result which 
satisfies the elastic unitarity relation or the phase theorem 0, is a lenghthy problem and 
will be dealt in a forthcoming publication [^T[. We would like to point out however the 
ambiguity of the corresponding integral equation, which is a generalisation of the MO 
equation, also exists but its nature is more complicated than that given in Eq. (|G|). 

We have shown in this note two methods to get agreement with experimental data. 
The first one consists in making a twice subtractions in the MO integral equation in order 
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to suppress the high energy contribution in the dispersion integral and hence it allow us 
to use the approximation of the elastic unitarity for the pion form factor. The second 
method is to take into account of the inelastic effect which, around the p region, can 
be simulated as the polynomial ambiguity of the solution of the MO equation. In both 
methods, the extra parameter is determined by the (low energy) rms radius of the pion 
and hence we can only calculate the p leptonic width and also its total width Eq. (It 
should be noticed that the rms pion radius can be calculated by a rapidly converging sum 
rule in terms of the magnitude of the pion form factor factor and the P-wave nir phase 
shifts using Eq. flD with n = 1 [p^]). 

Within the frame work of dispersion relation and unitarity which gives rise to the MO 
equation, there is no room for a constant contact term. Its existence could, however, be 
introduced as a large energy behavior of the form factor or some higher energy inelastic 
effect as explained in this note. 

This work is completed at the Institute for Nuclear and Particle Studies at the KEK 
High Energy Accelerator Research Organization. The author would like to thank Profes- 
sors M. Kobayashi and H. Sugawara for hospitality. 
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Figure Captions 

Fig.l : The phases of the function Q(s) or V(s) in degrees are given as a function of 
s(GeV 2 ). The experimental data are taken from references |0|, |16|, |PT|| . 

Fig. 2 :The square of the modulus of the function Q(s) is given as a function of s in 
GeV 2 (dashed line). The square of the modulus of the pion form factor V(s) as given by 
Eq. (|T3|) (solid line). Experimental data are taken from references []T3], [T4j. 
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